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Abstract— The variability of rainfall patterns has a significant impact on people's livelihoods especially in miren 
areas where rain-fed agriculture predominates. The temporal variability and trend of rainfall was analysed °° A oe 
using a newly rainfall time series created from the fusion with ground measurement and satellite data for s 
ten met stations for the period 1983-2021. Mann Kendall's (MK) non-parametric test was used to verify the thre aie 4 
rainfall trend, the variability was analysed using the Coefficient of Variation (CoV), the annual distribution ope 
of precipitation was checked using the precipitation concentration index (PCI) and finally the synchronicity : 
between stations was tested with an Kruskal Wallis (KW) test. The MK-test shows an increasing trend in QO) 
precipitation for all stations even if this trend is not significant. This increase resulted in a general increase 

in the amount of rain received in the area. The analysis of rainfall data revealed a moderate intra and inter 

annual variability (CoV+ 20%). The PCI also revealed an irregularity of rainfall due to its high seasonality. 

The analysis of variance highlighted an asynchrony of the stations. These results provide an in-depth 

understanding of the recent rainfall variation and trend in the cotton zone of Mali. 
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I. INTRODUCTION Located in West Africa (WA), Mali is a landlocked country 


Rainfall is one of the most key climate parameters which located populated with about 228 on of habitants 
regulate and determine farming activities and production (World Bank, 2022). The country is large with 1 245 000 


j : : : 
worldwide (Bekele et al. 2017; Mesike and Agbonaye, km/ and shares a border with seven countries (Algeria, 
Mauritania, Senegal, Guinea, Côte d’Ivoire, Burkina Faso 


and Niger). The country spans with four climatic zones from 
north to south (Sahara Desert, Sahelian, Sudanian and pre- 
Guinean). Agricultural activities including farming, herding 


2016). Agricultural boost in terms of high yield is based on 
sufficient amount and distribution of rainfall all through the 
seasons in any region (Zachariah et al. 2020). 
Consequently, in most African countries whose farming 
activity is mostly rainfed, robust and accurate estimation of 
the spatial and temporal distribution of rainfall including its 
trend are vital input parameters in order to secure 


and fishering, are the main activities and employ about 80% 
of the population and contribute about 36% to the GDP 
(Maiga et al. 2019). Given the fact that Malian agricultural 
activities are largely dependent on rainfed-agriculture, 
analysis of long-term historical rainfall data is imperative. 
These analyses will provide robust information on the 
distribution of precipitation on one hand and better prepare 


sustainable agricultural activities (Ayalew et al. 2012). The 
change in rainfall regime can’t be assessed easily due to 
meteorological processes but long-term rainfall could scale 
up the planning of agriculture in the rainfed region (Pradhan 
et al. 2020). 
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farmers on extreme events mitigation to avert agricultural 
losses on the second hand. 


The high inter-annual rainfall variability in space and time 
is one of the most relevant characteristics of West Africa in 
general and the southern zone of Mali in particular. This is 
exhibited by variable onsets of the rainy season, somewhat 
more predictable endings, and droughts or excess water 
occurrence at any time during the growing season. A 
number of studies on the contemporaneous debate (e.g. 
Traore et al. 2021; Sanogo et al. 2015; Giannini, 2015) 
show that rainfall increased in the Sahel following the 
severe drought of the 1980s. Conversely, other researchers 
believe that the Sahel experienced drought throughout the 
1990s (L’Hote et al. 2003; Nicholson et al. 2000). It is 
required to understand rainfall variability and trend at the 
temporal characteristics at local scale in order to gain a clear 
understanding of rainfall patterns. In this regard, long term 
study of rainfall pattern and its variability is imperative. 
Specially, in agricultural regions such as CPZ alteration in 
the rainfall patterns may have a major impact on the life and 
livelihood of people. With rapid climate change it becomes 
more imminent to study the trends in rainfall especially for 
those parts of the country which have higher dependability 
on agriculture. In this line, the present study seeks to 
investigate the recent variability and trend of rainfall and its 
distribution using 39-years (1983-2021) rainfall time series 
data for 10 selected weather stations in the cotton 
production zone (CPZ) of Mali. 


Il. MATERIAL AND METHODS 
2.1. Site description 


The cotton production zone (CPZ) of Mali is located 
between 3°59'02" and 10°36'09" West longitude and 
10°09'45" and 14°23'13" North latitude (Figure 1) with a 
coverage of 150,000 km? and an estimated population of 8 
million habitants in 2017 (Soumare and Traore, 2019). The 
CZ extends from the Sahelian (300-700 mm) in the north, 
Sudanian (700-1200 mm) in the centre to the Sudano- 
Guinean (1200-1600 mm) eco-climatic zone in the south. 
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The spatial distribution of vegetation is largely related to 
cumulative (annual) rainfall and the length of the rainy 
season, which varies along the eco-climatic gradient. 
Vegetation patterns on a landscape scale are determined by 
localised climatic variations and human activities such as 
bush clearing and deforestation for agricultural or energy 
purposes, overgrazing and gold mining. The cotton 
production in the area is led by the Malian Company for the 
Development of Textiles (CMDT), spatially organised into 
4 subsidiaries, 41 sectors and 520 zones of agricultural 
production involving more than 4,000 villages. Smallholder 
farming is the most common agricultural practice in the 
region. Millet, sorghum and cotton are primarily cultivated, 
whilst pastoralism (bovines, goats and sheep) is practised 
throughout the area. Cotton occupies about 30% of the 
cropping area, 60% for cereals (maize, sorghum and millet), 
and only 10% of the land is given over to small crops 
(cowpeas, soybeans, groundnuts, etc.). 


2.2. Data description 


The World Meteorological Organization (WMO, 2018) 
recommends an observation period of more than 30 years to 
ensure the independence of climate data time series to cope 
with natural climate variability. This fact is because shorter 
time series are more sensitive to the values at the start and 
end of the series. Therefore time-series data of rainfall from 
1983 to 2021 for ten (10) meteorological stations in the 
cotton production zone were created. These data were 
derived from a fusion of satellite rainfall estimates from 
African Rainfall Climatology version 2 (ARC-2) adjusted 
with weathers stations data using the period of overlap of 
the two data sources then extended over the study period. 
Calibration procedures are described in Traore, ND.; Traore 
et al. (2022). Table 1 presents the general geographic 
information of the selected weather stations. These stations 
called reference weather stations in this study are those 
among the 54 stations listed in the area with no more than 
two (04) years of missing data. This is the criterion that 
motivated their choice for the present work. 
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Fig.1: Relief map of cotton zone with selected weather station’s location within the cotton zone of Mali (Source: Traore et al. 
2022) 


Table 1: Selected weather stations of the cotton zone and their general geographic information. 


Long CW) Lat (N) Altitude 
Code Station name AEZ* DDMMSS DDMMSS (m) 
270131 Senou PK 7°57'00" 12°31'59" 373 
270162 Bougouni HBN 7°28'35" 11°25'34" 350 
270141 Dioila PK 6°46'32" 12°29'14" 309 
270178 Kadiolo HBN 5°45'50" 10°33'34" 353 
270107 Kita PM 9°29'43" 13°02'34" 349 
270144 Koutiala PK 5°27'38" 12°18'35" 411 
270121 N'tarla IRCT PK 5°45'00" 12°42'00" 327 
270100 San BH 4°54'03" 13°18'03" 284 
270165 Sikasso PK 5°39'26" 11°19'25" 355 
270130 Sotuba PM 7°55'589" 12°39'00" 315 


AEZ: Agro-ecological zone; PK: Plateau of Koutiala, HBN: Haut Bani Niger, PM: Mandingo Plateau, BH: Plateau of 
Bandiagara 


ISSN: 2456-1878 (Int. J. Environ. Agric. Biotech.) 
https://dx.doi.org/10.22161/ijeab.86.5 38 


Traore International Journal of Environment, Agriculture and Biotechnology, 8(6)-2023 


2.3. Analyse procedures 
Normality test 


Assessment of the normality of data is a prerequisite for 
selecting the best statistical methods for data analysis. In 
this regard, the available two main methods for assessing 
normality are graphical and numerical including statistical 
tests (Machin et al. 2007; Bland et al. 2015). This study 
applied a normality test using the Kolmogorov-Smirnov test 
(KS-test) method. KS-test is recommended when sample 
size is greater than 50 (Mishra et al. 2019). In such, daily, 
monthly and annual rainfall for all the stations were 
subjected to normality tests. The null hypothesis states that 
data are taken from a normal distributed population and 
when P > 0.05, null hypothesis is accepted and data are 
called as normally distributed. 


Coefficient of Variation (CoV) 


The study uses CoV to examine the variability in rainfall. 
The higher value of CoV is an indicator of larger rainfall 
variability, and vice versa. The equation of CoV is as 
follows: 


CoV = 2x100 (1) 


Where CoV is the coefficient of variation; O is the standard 
deviation and Æ is the mean of rainfall over the studied 
period. CoV was computed to classify the degree of rainfall 
variability according to El-Mahdy, 2021 and Hael, 2021. 
CoV less than 20 indicate low variability, CoV greater than 
20 and less than 30 indicate moderate variability; and CoV 
greater than 30 indicate high variability (Royé and Martin- 
Vide, 2017). 


Mann Kendall trend test 


From the wide types of trend analysis methods, this study 
uses the non-parametric Mann-Kendall trend test (Mk-test) 
proposed by Mann (1945) and Kendall (1975). This test 
does not require the data to be normally distributed and has 
low sensitivity in abrupt breaks due to inhomogeneous 
time-series (Tabari et al. 2011) and have been widely 
employed to detect monotonic trends in the time series of 
hydrometeorological variables (Asfaw et al. 2018; Xu et al. 
2018; Zakwan and Ara, 2018; Mohamed and El-Mahdy, 
2021; Mohamed et al. 2022). The MK test is grounded on a 
null hypothesis (Ho), which indicates that there is no trend- 
the data are independent and randomly ordered-and this is 
verified against the alternative hypothesis (Ha), which 
supposes that there is a trend (Koudahe et al. 2018). In the 
calculation, the value of Z is the judgement criterion for the 
trend change (Xu et al. 2018). When |Z| < 1.96, the null 
hypothesis HO is accepted, indicating that there is no 
significant trend at the 0.05 significance level. |Z| > 1.96 
demonstrates the trend of the time series is statistically 
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significant. It must be noted that a positive Z indicates that 
the sequence has an increasing trend, while a negative Z 
reflects a declining trend. 


Precipitation concentration index 


Precipitation Concentration Index (PCI) was used to assess 
the monthly heterogeneity of rainfall amounts. PCI is a 
useful indicator to determine the precipitation changes of a 
specific region and defined as the ratio between sum of 
squared monthly rainfall to the square of annual rainfall 
(Bhattacharyya and Sreekesh, 2022). The PCI, proposed by 
Oliver (1980) and further modified by De Luis et al. (2011), 
is used for the calculation of the annual PCI following: 
Zii Pi 


PCI = 100 x SG 


(2) 


PCI values of less than 10 indicate quite uniform annual 
distribution of rainfall, values between 11 to 15 denote a 
moderate rainfall distribution, 16 to 20 denote irregular 
rainfall distribution and above 20 represent a strong 
irregularity of rainfall distribution Asfaw et al. 2018 and 
Rahman et al. 2019, Pawar et al. 2022. 


Comparative analysis between stations 


Comparison of several groups or independent samples 
requires the application of an appropriate statistical test. In 
this regard, the common statistical tests are ANOVA for 
normally distributed samples or groups and the Kruskal 
Wallis test (KW-test) for non-normally distributed samples 
or groups (Cabral Júnior and Lucena, 2019). In this study, 
the KW-test test was used since the hypothesis of the 
normality distribution of data was rejected at 1% of 
Statistical significance, as verified by the KS-test (p-value: 
0.5). As the KW-test compares (paired or unpaired) k 
samples based on the null hypothesis that the median 
differences within groups are not significant. In this 
research the groups were formed by daily/annual rainfall 
data for individual met stations to check the rainfall pattern’ 
similarity and difference (synchronicity and asynchronicity) 
during the 39 years period. The null hypothesis is that there 
are no significant differences between the rainfall medians 
in the pattern in the study area. 


Il. RESULTS AND DISCUSSION 
3.1. Long-term rainfall temporal variability 


Statistics of annual rainfall (1983-2021) show a 
comparatively higher rainfall amount for all the stations. 
The mean annual amount of rainfall varies from 1223.3 mm 
with a standard deviation of 272.6 at Sikasso to 745.4 mm 
with a standard deviation of 155 at San. The maximum 
rainfall at Sikasso was 1848.3 mm recorded in 2018 and the 
minimum was 781.0 mm recorded in 2002. In San, the 
maximum rainfall of the period was 1135.4 mm recorded in 
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2020 and the minimum was 480.9 mm recorded in 1992. 
The result of the normality test using KS revealed non- 
normality in the series for all the ten stations at a 
significance of 5%. Rainfall amounts show low variability 
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for Kadiolo station (18.8%) and Bougouni station (19.4%). 
Moderate variability in rainfall is observed in the rest of the 
stations. 


Table 2: Summary statistics of annual rainfall (mm) of selected weather stations (1953-2021) 


Station Min Median Max Average Standard deviation CoV (%) 
Sikasso 781.0 1155.4 1848.3 1223.3 258.7 21.6 
Bougouni 807.5 1151.7 1790.0 1187.7 227.8 19.4 
Kadiolo 657.4 1115.5 1892.5 1192.5 218.7 18.8 
Dioila 539.9 1024.7 1796.8 1074.1 271.3 25.7 
Senou 577.9 986.5 1597.8 1017.7 207.1 20.4 
Kita 630.2 978.9 1726.7 1033.4 226.2 22.7 
N’tarla 557.4 957.5 1608.0 977.4 236.1 24.1 
Sotuba 535.6 941.3 1562.3 982.9 199.5 20.5 
Koutiala 588.8 912.0 1587.1 953.5 261.9 26.6 
San 480,9 719.0 1135.4 745.4 153.4 20.7 


3.2. Rainfall trend 


The MK test and Sen’s slope had been calculated for the 
trend analyses on the long period of rainfall data from 1983 
to 2021 for the 10 meteorological stations in the zone. The 
results of trend analyses by the MK-test for the study 
stations were displayed in Table 2. In the Mann Kendall 
trend test, a p-value less than 0.05 means a significant trend 
and a p-value greater than 0.05 means is considered an 
insignificant or simply no trend (Di Leo et al. 2020; 
Amrhein et al. 2019). In general, all the stations showed a 
positive Sen’s Slope (between S Slope= 3.99 and 19.67) 
indicating at least an increasing rainfall trend. Kita in 
western zone is the only station which exhibits a non- 
statistically significant trend (p-value of 0.345 greater than 


0.05) with a positive Sen’s slope value (Sen’s slope = 3.99). 
The other stations show p-values between p=0.000 and 
p=0.021 (p < 0.05) which makes the trends statistically 
significant upward trend of total rainfall amount. All the 
stations show either zero or positive values for Kendall’s 
tau, Mann Kendall parameter (S), Variance (S) and this is a 
good sign of the increasing trend of rainfall. These results 
agreed with the rainfall recovery trend asserted by Lalou et 
al. 2019; Giannini, 2015; Sanogo et al. 2015 in West Africa. 
Indeed, the results are also in agreement with the analysis 
of Traore et al. 2021 who had noticed an increased rainfall 
trend for the period 1983-2018 using a rainfall time series 
from one meteorological station in the CPZ of Mali. With 
these results, one can assume that there is an increasing 
trend of rainfall amounts over the CPZ, significant at 99%. 


Table 3: Mann-Kendal rainfall trend and Sen Slope for selected weather station (1983-2021) 


Kendall's Sen Slope Z 

Station tau S Var(S) P-value Alpha HO Value 
Sikasso 0.34 12.01 250 6832.7 0.003 0.05 Rejected 3.01*** 
Bougouni 0.26 8.45 192 6832.7 0.021 0.05 Rejected 2.31** 
Kadiolo 0.32 10.44 238 6832.7 0.004 0.05 Rejected = 2.87*** 
Dioila 0.58 19.67 432 6832.7 < 0.0001 0.05 Rejected  5.21*** 
Senou 0.46 12.04 338 6832.7 < 0.0001 0.05 Rejected 4.08*** 
Kita 0.11 3.99 79 6831.7 0.345 0.05 Accepted 0.94 

N'tarla 0.57 15.44 420 6832.7 < 0.0001 0.05 Rejected  5.07*** 
Sotuba 0.34 8.43 250 6832.7 0.003 0.05 Rejected 3.01*** 
Koutiala 0.41 11.83 302 6832.7 0 0.05 Rejected 3.64*** 
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San 0.43 8.58 320 6832.7 0 0.05 Rejected 3.86*** 


Significance level: *** significant at 0.001; ** significant at 0.05 


3.3. Precipitation Concentration Index (PCI) Analysis 


The variation and distribution in seasonal rainfall signify a 
significant indicator to appreciate the concentration and 
sequential distribution of precipitation over a given year. 
Therefore, the variability and concentration of rainfall for 
seasonal rain was evaluated using annual PCI based on the 
monthly rainfall over a 39-year period (1983—2021) for 10 
meteorological stations in the cotton zone of Mali. The 


result shows that the rainfall at the 10 weathers stations falls 
under two classes of moderate and irregular precipitation 
concentration. Kadiolo station with a PCI value of 15 is the 
only station highlighting moderate rainfall distribution. A 
PCI value between 16 and 20 were observed in the reminder 
stations denoting an irregular rainfall distribution. These 
results corroborate with the results of Quenum et al. 2020 
who had already proven an irregularity of rainfall in the 
savannah zone because of its more pronounced seasonality. 


Table 4: PCI distribution for the ten stations used during 1983-2021 period 


PCI distribution PCI range Annual 
Uniform precipitation < 10 NA 
Moderate precipitation 11 to 15 Kadiolo 
Sikasso, Bougouni, Dioila, Senou, Kita, N’Tarla 
I 1 1 t 2 3 2 3 2 td d 
ee eta! Sotuba, Koutiala, San 
Strong irregularity > 20 NA 
2000 
A A A 
AB AB 
1750 
B B B 
c 
1500 
- 1250 
a 
750 p 
500 
Sikasso Bougouni Kadiolo Dioila Kita Ntarla Sotuba Koutiala San 


Fig.2: Boxplot classification of groups weather stations 


3.4. Inter-station comparison of rainfall 


Figure (3) shows the result of the comparison of the rainfall 
series of the 10 weather stations. Three homogeneous 
groups of stations were obtained whose medians are 
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different from each other. The group “A” consists of the 
stations of Sikasso, Bougouni and Kadiolo with a long-term 
average rainfall of around 1200 mm with a standard 
deviation of 253. The group “B” consists of the stations of 
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N’tarla, Sotuba and Koutiala with a long-term average 
rainfall of around 971mm with a standard deviation of 231. 
The “C” group consists only of the station of the station of 
San with a long-term average rainfall of 745 mm with a 
standard deviation of 155. The stations of Dioila, Senou and 
Kita, intermediate group “AB", present almost the same 
characteristics of the group “A” and the group “B” with a 
long-term average rainfall of about 1042 mm with a 
standard deviation of 252. Through this result, one can 
understand that there is local synchronicity between the 
stations and not zonal. This means that an increase or 
decrease in rainfall in any weather station is not widespread 
for all reminder weather stations in the same climatic zone 
or even all stations in the study area. 


IV. CONCLUSION 


In this paper, an investigation of the rainfall variability and 
trend in the CPZ of Mali has been performed by the means 
of daily to annual rainfall dataset from 1983 to 2021. 
Statistical analysis including CoV, Mann-Kendall non 
parametric trend test and PCI were carried out in order to 
detect variability, possible trend and concentration in the 
rainfall time series data. The rainfall amount showed 
moderate annual variability. A significant positive trend in 
total rainfall was observed over the study time period. The 
PCI showed an irregular rainfall variability in general that 
corresponds to the characteristic of the savannah zone. The 
comparison analysis used highlights a synchronicity 
between some weather stations. Finally, for improving 
precision and reliability of the application of the findings 
for practical use, increasing the number of study weather 
stations would be crucial. However, this research shows the 
importance of local level study in practical decision-making 
processes in agriculture and other water management 
programmes in the climate change adaptation context. 
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